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1. Introduction 

In ultra-relativistic heavy-ion collisions [1, 2] large numbers of pions can be created 
in small space-time regions. In such a situation we may ask the question: What is the 
probability that a part of them will form a condensate? In the recent paper of Greiner, 
Gong and Mtiller [3], it has been pointed out that the occurrence of condensed pions can 
only take place if the system evolves far away from thermal equilibrium. They argue that if 
the hadronization of a deconfined quark-gluon plasma (QGP) occurs rather slowly and the 
initial pion temperature is of the order of 170 MeV, then the inelastic reactions changing 
the pion number enforce the pion chemical potential to be zero. The detailed study [3] 
of the subsequent adiabatic expansion of such a pion gas indicates that the pion chemical 
potential builds up. However, it remains always below the critical value, independent of the 
interactions assumed. Consequently, the pion gas in thermal equilibrium can never reach 
the onset of Bose- Einstein condensation (BEC). 

Having the above results in mind, in the present paper we are going to consider a non- 
equilibrium hadronization of the QGP and to study whether already during such a process a 
part of pions can produce a condensate. To some extent, our investigations can be regarded 
as supplementary to those of [3] since we concentrate on possible pion condensation occurring 
at earlier times and employ a different approach, i.e., we solve the system of kinetic equations 
in a model study, that allows us to study the non-equilibrium processes. 

In our calculation, we assume that the hadronization is fast, so, at least for the initial 
stages of the evolution of the system, we can assume that it happens within a constant 
volume. Nevertheless, the produced pions can leave this volume (evaporate), which is crucial 
for the overall increase of the entropy during this reaction. One can check that the total 
hadronization of the QGP into the pion gas taking place in a constant volume and at constant 
energy is impossible, since, roughly speaking, it leads to the reduction of the number of 
degrees of freedom (if one takes into account only pions and neglects resonances). On the 
other hand a partial hadronization is possible, and, if connected with the simultaneous 
evaporation of mesons, can lead in a sequence of steps to the complete hadronization of the 
plasma and to the increase of entropy. 

We also assume that the reactions conserving the pion number dominate, i.e., except for 
the very process of the hadronization, the other reactions do not change the pion number. 
Such an approach can be justified if the cross sections for the inelastic, pion number changing 
reactions such as tttttt tttt are much smaller than those for elastic ones. Since the in- 
medium cross sections are rather poorly known, we think that it is useful to look for the 
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consequences of this hypothesis. In such a case, a very dense system of pions is created 
just at the beginning of the hadronization process, and, perhaps, a part of them will form a 
condensate. 

In order to study quantitatively the possibility of appearance of the pion condensate, we 
shall solve a system of kinetic equations describing the transformation of quarks and gluons 
into pions. To facilitate these considerations, we assume that the phase space distribution 
functions are functions of momentum p = |p | and time t only. Thus the particles are always 
uniformly distributed in space and the system is isotropic. Moreover, we consider the kinetic 
equations in the relaxation time approximation, reducing the problem to solving a system 
of intego-differential equations. 

A novel feature of our approach is the incorporation of the phenomenon of the BEC into 
such a relaxation time formalism. To some extent, our description is close to the statistical 
one: For a given pion distribution function i), we check whether the corresponding 

thermal distribution (i.e., the distribution which has the thermal shape and gives the same 
energy and particle densities as the function /^(p, i)) contains a condensate part. If so, then 
the pion collisions that always bring the distribution function /^(p, i) to the thermal one 
with the speed characterized by the relaxation time, will try to force some part of pions to 
form the condensate. 

Due to the pion interactions between themselves and with quarks and gluons, the pion 
dispersion relation in medium can change. In order to take into account some aspects of 
this problem, we shall consider different pion masses in our calculations, i.e., we shall take 
the pion dispersion relation of the (free) form = + with different values of m^. In 
particular, we shall consider the cases = and = 140 MeV. 

The paper is organized as follows: In the next Section, we introduce the kinetic equations 
in the relaxation time approximation, define the thermal functions, discuss the increase of 
entropy, and calculate the relaxation times. Section 3 contains our results. We summarize 
and conclude in Section 4. 
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2. Kinetic description of the hadronization of the QGP 



i) Kinetic equations in the relaxation time approximation 

We assume that the hadronization reactions take place in constant volume Vpi and that 
the system is spherically symmetric, i.e., Vpi = ^T^Rpi where Rpi is the radius of the fireball 
containing quarks, gluons and pions (quarks and gluons cannot escape from this mixed phase 
region). The number of pions evaporated per unit time can be estimated to be \f-K{jPit)SpiVp, 



where 5"^; is the area of the fireball surface and Vp = 2p/Tvyp^ + average pion 

velocity perpendicular to the surface. The factor | accounts for the fact that only one half 
of the particles situated close to the surface and having momentum p = |p | move outwards 
and can, in consequence, leave the system. The volume occupied by the emitted pions is 
defined as Vhad{t) = f^r [-R^ad(i) — -Rp/] where Rhadif) is the radius of our total hadronic 
system. We assume that Rhadif) increases in time with the velocity of light. 

Having this geometry in mind we study the following set of the kinetic equations (they 
are a generalization of the equations introduced in [4] by including gluons and, as we shall 
see later, by considering the quantum statistics) 



dfq{P,t) ^ fq{P,t) - fq,th{P,t) _ fq{p,t) 
dt Tth{t) Tq^had{t) ' 

dfg{p,t) ^ fg{p,t) - fg,th{p,t) _ fg{p,t) 
dt Tth{t) Tg^had{t) ' 



(1) 

(2) 



dfAP^ ^ fAP^t) - U,th{p,t) ^ fq{p,t) ^ ^ fg{p,t) ^ _ y 

dt Tth{t) Tq^had{t) Tg^had{t) 2Vpl 

dfem{p,t) fem{p,t) dVhad{t) . Spl _ 

dt - VHad{t) dt ^ ^""^^' ^ 2VHad{tf''- ^ ^ 

Eq. (1) determines the time evolution of the quark distribution function fq{p,t) (it 
describes both quarks and antiquarks, i.e., fq{p,t) = fquarks{p,t) + fantiquarks{p,t)). The first 
term on the right-hand-side (RHS) of Eq. (1) is the collision term written in the relaxation 
time approximation; it is responsible for the thermalization of quarks since the distribution 
function fq{p,t) is always attracted to the thermal one fq,th{jp-,i)- The second term describes 
the loss of quarks due to the hadronization process. We note that the number of quarks in 
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the system will be always equal to the number of antiquarks. Consequently, we study the 
hadronization of the baryon-free plasma and this might lead to the creation of the condensate 
which carries no isospin. 

Eq. (2) has the form completely analogous to that of Eq. (1) and describes the ther- 
malization and hadronization of gluons. In Eq. (3), for the time evolution of the distri- 
bution functions of pions, the first three terms on the RHS describe pion thermalization 
and pion creation due to the hadronization of quarks and gluons, respectively. The factors 
Rq = + ml)/{p^ + ml) and Rg = ^(p^ + ml)/{p^ + ml) (here m^^mg and is the 

quark, gluon and the pion mass) are included in order to have the overall conservation of 
the energy in the system. The fourth term accounts for the evaporation of pions. 

Our last kinetic equation (4) describes the time changes of the distribution function of 
the pions which have been emitted from the reaction volume. The first term on the RHS of 
Eq. (4) appears due to the expansion of these pions into the vacuum, and the second term 
collects the pions acting as a sink in Eq. (3). 

The effects of Bose- Einstein condensation can be included into our description if we 
represent the pion distribution function as a sum of two terms, namely 

MP,t) = f>{p,t) + c.(i)(27r)V^)(p). (5) 

The analogous decomposition holds also for the function fT^^th{jP-,i)- In this case Eq. (3) can 
be rewritten as follows 



df^M f^M-fl,HM fg{p,t) fg{p,t) f>.,.s,i 

at Tth{t) Tq^had{t) Tg,had\l') ^Vpl 

dCT,{t) _ CT,{t) - CT,^th{t) /-X 

dt ~ nnit) ' ^ ^ 

where we have explicitely separated the smooth and the singular parts of the pion distribu- 
tions. 

ii) Thermal distributions 

Eqs. (1) - (4) can be solved for any initial conditions (i.e., assuming some particular 
form of the distribution functions at the initial time t = 0) provided the thermal distribution 
functions fq,th{jp-,t)-, fg,th{jp-,t) and f-K,th{jP-,t) are known at all times. Since the first terms 
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on the RHS of Eqs. (1), (2) and (3) represent the collision terms (in the relaxation time 
approximation) they must obey the symmetry leading to the energy conservation. This gives 
the following constraint for fq,th{jP->^)-> fg,th{jP->^) f-K,th{jP->^) each time t 



where i = q,g ot tt. Our collision terms for quarks and gluons account for elastic and inelastic 
scattering, which leads to the change in the number of participating quarks and gluons. On 
the other hand, in the case of pions, we assume that their number is not changed by collisions 
[except for the hadronization processes, but these are already included in the separate terms). 
Consequently, the pion collision term should posess the additional symmetry, namely 



/ 



_fp_ 
(2^ 



(9) 



Eqs. (8) and (9) can be used to determine the thermodynamic parameters appearing in 
the thermal distributions. For quarks and gluons the latter have the form 



and 



fq,th{P^ 9ql exp 



fg,th{P^ = 9g\ exp 



p2 + 



P 



+ 1 



(10) 



(11) 



where we have introduced the degeneracy factors: gq = 24 (quarks and antiquarks having 
two different spin projections, 3 colors and 2 flavors) and gg = 16 (two spin projections times 
8 colors). In the case of pions the thermal background distribution function can be written 
as 



-l| +c.,t,(i)(27r)¥^)(p), (12) 

where 5^ = 3 (three different values of the isospin). For pions we have generalized the 
thermal distribution in order to include also the effect of the condensation. The number of 
pions in the thermal condensate equals c^,fh(i), and //^(O pion chemical potential. The 



f-K,th{P^ = 9-k{ exp 



p^ +ml- fi^{t) 



6 



quantities Ti[t) appearing in Eqs. (10) - (12) are the quark, gluon and pion temperatures, 
respectively. 

In the situation without the pion condensation Eqs. (8) and (9) can be fulfilled by finding 
the appropriate values of Ti[t) and /^^(i), and by setting c^(i) = c^,fh(i) = 0. For increasing 
pion densities the pion chemical potential //^(O increases and at some stage (let us say for 
t — tcond^ it may happen that it reaches its maximal value, i.e., it becomes equal to the 
pion effective mass m^. This is exactly the moment when the pion condensate sets in. For 
larger pion densities we have to assume that a part of them occupies macroscopically the 
zero momentum state. In this situation Eqs. (8) and (9), considered for i = tt, are two 
equations determining the pion temperature Tn-(i) and the pion thermal condensate c^,fh(i) 
(the pion chemical potential for t > tcond becomes constant and takes on its critical value 
/i^ = m^). If c^,fh(i) starts out as being different from zero, then according to our Eq. (7), 
the actual distribution function /^(p, i) also develops a singular part c^(i). 

Hi) Energy conservation law and the increase of entropy 

Our kinetic equations (1) - (4) are constructed in such a way that the total energy E[t) 
of quarks, gluons and pions (in and out of the fireball) is conserved. Defining the particle 
energy densities as 

ei{t) = j j^Jp'+mUi{p,t), (13) 
we can explicitly write this conservation law in the form 



dt 



^(^) = Ji{^pi £,{t) + eg{t) + e^{t) +Vhad{t)eem{t)^ = 0. (14) 



The total entropy of our mixture can be defined as follows 



S{t) = Vpl Sg{t) + Sg{t) + S^{t) 



+ Vhad{t)Sem{t), 



(15) 



where [5] 



J (27r)3 



d'p ifi{p,t)^^fi{p,t) 1 



1 + e, 



fi{p,t) 



In 



1 + e, 



fi{p,t) 



(16) 



Here the quantity equals 1 for gluons and pions, and —1 for quarks (the classical Boltzmann 
definition corresponds to the limit e ^ 0). It is important to notice that in the case of pions 
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the distribution function /^(p, i) in (16) should correspond to f^[p,t) since the condensate 
does not contribute to the entropy of the system. 

Although Eqs. (1) - (4) lead to the conservation law (14) they do not guarantee, in 
general, that their solutions correspond to states of increasing entropy. Consequently, for 
each solution we check the behaviour of the function S{t), and accept as the physical ones 
only these for which dS{t)/dt > (for each t). This fact indicates some difRciency of the 
present description of hadronization, namely it cannot be valid for arbitrary initial conditions 
or parameters. One might expect this type of difficulty, since some physical processes have 
been neglected in our approach and they might become important for a given class of initial 
conditions. In particular, we have neglected processes such as: (i) an expansion of the volume 
of hadronizing matter, (ii) the possible back reactions of pions into quarks and gluons, (iii) 
the production of pion resonances, or (iv) the process in which a quark tries to leave the 
plasma, stretches the string which later decays producing the meson outside. One observes 
that all these processes decrease chance of producing of a condensate. 

The study of the time evolution of the entropy shows, in particular, that it is crucial to 
take into account the appropriate balance between the hadronization and evaporation rates, 
i.e., the hadronization cannot be too fast with respect to evaporation. In our case, we control 
the speed of hadronization by considering different values of the hadronization cross section 
cr, which is regarded as a parameter. The evaporation rate is determined mainly by the 
geometry of the system, in particular, by the surface to volume ratio. As a consequence, we 
can always study opizma/ conditions for creation of the condensate (i.e., for a given geometry 
we consider the largest value of a which allows for entropy growth) and check whether it has 
been really produced. 

iv) Relaxation times and cross sections 

Let us now turn to the discussion of the relaxation times appearing in our kinetic equa- 
tions (1) - (4). The analysis of the exact collision term in the Boltzmann kinetic equation 
(analogous to that in [6]) leads to an expression for the average time for hadronization of 
two quarks into two pions 



where crqg_>^^(5) denotes the total cross section for this process. The relativistic flux factor 




(^) 



(17) 



rg,had{P,t) 2^p2 ^^2 J (27r)3y5[T^ 



of incoming quarks is Fqq(s) = ^^s(s — 4mg), with ^/s being the center-of- 



"-mass energy of 
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the quarks with momenta p and pi, respectively. The analogous expression holds, of course, 
for the gluon hadronization time Tg^had{jP->^)- 

Here we also make simplifying assumptions: In the limit of vanishing masses (this case is 
exact for gluons and represents a very good approximation for current quarks) and for the 
constant cross sections (cr's independent of s) we find 

1 _ 1 _ 1 _ 1 

Tq,had{t) ^ Tg^had{t) ^ 

where 

ni{t) = j-^JiM (19) 

are particle's number densities. Consequently, one finds that the hadronization relaxation 
times are functions of time only. We note that if the cross sections depend on 5, the relaxation 
times (even for vrii = 0) can exhibit a p-dependence. In particular, the singularities of the 
cross sections can be reflected in the singularities of the inverse relaxation times. For the 
discussion of this point in the context of the critical scattering we refer the reader to Ref. 
[4]. 

The relaxation time for the thermalization of quarks (determined mainly by the elastic 
collisions) can be written in the form 



/ ^7^^^=/7r(pi, 0^9^(5)0-5^^5^(5), (20) 



where the three terms on the RHS of Eq. (20) appear due to the quark-quark, quark-gluon 
and quark-pion collisions, respectively. The formulae for the thermalization times of gluons 
and pions can be obtained from Eq. (20) by appropriate changes of the indices q,g and tt. 
Considering again the massless limit and assuming that the cross sections are independent 



of 5, we find 



7: [nq{t)aqq^qq + ng{t)(7qg^qg + n^{t)(7g^^g^] . (21) 



rth{t) 2 

In order to reduce the number of independent parameters we shall set all the cross sections 
appearing in Eqs. (17) and (20) to be the same and equal a. In this case the thermalization 
time is the same for quarks, gluons and pions. 
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We note that in our calculations we shall consider different effective pion masses, in 
particular, we shall take into account the normal vacuum mass = 140 MeV. In such a 
case the approximation is not as good as the approximation 0. Nevertheless, 

we have checked that the inclusion of the finite meson mass here (i.e., in the calculation of 
the relaxation times) represents only a small correction. Much more important are effects 
coming from the appearance of the finite and large (compatible with the temperature) pion 
mass in the thermal functions (12) since in this case is the argument of the exponent; 
these corrections will be exactly taken into account. 

3. Results 

The initial conditions assumed for solving the kinetic equations (1) - (4) correspond to 
a quark-gluon plasma in thermal and chemical equilibrium, which is confined in a spherical 
fireball. This means that the quark and gluon distribution functions have initially the ther- 
mal form (10) and (11). Moreover, at the beginning of the process there are no pions in the 
system, i.e., U{jp,t = 0) = 0. 

For i = we specify the initial plasma temperature To and the plasma radius Rpi (initially 
yhad = 0). The additional parameters are: the value of the hadronization/thermalization 
cross section cr, and the quark and pion masses (the gluon mass is always taken to be zero). 
We do the calculations for different values of the pion mass since we want to study the 
sensitivity of our results to the change of the effective pion mass resulting from the pion 
interactions with other particles in the medium. 

Let us first consider the case when all the masses are equal to zero [vriq = rUg = = 0). 
In Figs. 2 and 3 we show our results for the initial temperature To = 150 MeV, a = 9 mb 
and Rpi = 2 fm [Vpi ~ 34 fm^). This situation corresponds to the initial energy density of 
0.8 GeV/fm^, and the total number of particles (quarks and gluons) is 59. The initial total 
entropy calculated from (15) is 227. The choice To = 150 MeV is motivated by the results of 
the lattice simulations of QCD, which indicate such a temperature of the phase transition. 

In Fig. 1 we plot the time evolution of the thermodynamic parameters characterizing the 
pion thermal distribution. The upper curve describes the changes of the temperature and the 
lower one shows the behaviour of the chemical potential. The pion temperature first increases 
and later saturates at 200 MeV. In our opinion, the increase of the pion temperature is due 
to the fact that we do not have pion number changing processes. (At the same time the 
quark and gluon temperatures remain approximately constant.) The pion chemical potential 
is initially very large and negative since the pion density is very low, however later it builds 
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up very fast and for t = 0.4 fm it already reaches its critical value, i.e., in this case zero. 

In Fig. 2 the time evolution of the particle's numbers has been drawn. We can see that 
the numbers of quarks and gluons decrease in time. Initially the number of pions increases 
rather fast since the process of evaporation (proportional to /^(p, i)) is rather weak. At the 
moment when the pion chemical potential takes on its critical value, some fraction of them 
starts to occupy the zero momentum state, i.e., for t > 0.4 fm we find Ncond ^ 0. For larger 
times the number of pions in the fireball saturates, however the number of evaporated pions 
starts to increase rapidly (at this stage the evaporation becomes a very efficient process). 
After 2 fm the average number of pions in the condensate is 2. 

In the considered case (corresponding to Figs. 1 and 2) the entropy increases only by 
7%. In fact, if we did the calculations for larger values of the hadronization cross sections, 
we would find that the entropy derivative becomes negative in some time intervals. In this 
place we are touching the problems already discussed before: For some initial conditions 
(or values of the parameters) one should include more processes into the kinetic descrip- 
tion. Nevertheless, we expect that such improvements can have only little effect on general 
conditions concerning the condensate creation. Our specific calculation shows that for very 
large a (larger chances for condensate cration) we have to take into account other processes 
(reducing the probability of pion condensation). Consequently, the overall result of more 
complicated calculations can be very similar to the present ones. 

Let us now discuss the results obtained for the following set of masses: rUq = 10 MeV, 
rUg = 0, and = 140 MeV. The initial temperature is again equal to 150 MeV and Rpi = 
2 fm. The hadronization cross section is 7 mb since for larger values the total entropy would 
be (in some time intervals) a decreasing function of time. One can notice that in this case 
(massive pions) the hadronization process cannot be so fast as in the previous case since 
the evaporation is slower. In Fig. 3 (analogous to Fig. 1) we show the time evolution of 
the thermodynamic parameters: The pion temperature slightly increases, whereas the pion 
chemical potential grows up substantially. Nevertheless, its value saturates for larger times 
and remains always below the critical one (in this case equal to 140 MeV). Consequently, in 
this case pions do not condense. 

4. Conclusions 

In this paper we have studied the hadronization of a QGP into the pion gas and inves- 
tigated the possibility of creation of the pion condensate during such a process. For fixed 
initial temperature we find that the results depend on the pion effective mass in medium. 
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For massless pions, a small fraction of them will condense. On the other hand for pion 
effective masses compatible with its vacuum mass, no onset of Bose- Einstein condensation 
is found. The requirement that the entropy should increase during the hadronization leads 
to the strong constraints on the possibility of creation of a very dense pion system and, in 
consequence, on possibility of pion condensation. 

We studied the initial stages of hadronization since for larger times other processes should 
become important (e.g., pion number changing reactions or volume expansion) which has 
been neglected in our approach. Such processes will diminish the chances of producing a 
condensate and, consequently, if it is not produced fast then it is probably not produced at 
all. Another open question is whether a condensate once produced has chances to survive. 
In our case we observe the formation of the condensate for massless pions. The effect is not 
large and, perhaps, can be destroyed when the particles gain their true vacuum masses. 
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Figure Caption 

Fig. 1 Pion temperature and pion chemical potential as the functions of time (To = 150 
MeV, a = 9 mb, rUq = rUg = = 0, and Rpi = 2 fm). 

Fig. 2 Time evolution of the particle's numbers. The parameters as in Fig. 1. 

Fig. 3 Pion temperature and pion chemical potential as the functions of time (To = 150 
MeV, a = 7 mb, = 10 MeV, = 0,m^ = 140 MeV, and Rpi = 2 fm). 
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